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We  study  the  temperature  dependent  Young’s  modulus  for  the  glass/ceramic  seal  material  used  in  solid 
oxide  fuel  cells  (SOFCs).  With  longer  heat  treatment  or  aging  time  during  operation,  further  devitrification 
may  reduce  the  residual  glass  content  in  the  seal  material  while  boosting  the  ceramic  crystalline  content. 
In  the  meantime,  micro-voids  induced  by  the  cooling  process  from  the  high  operating  temperature  to  room 
temperature  can  potentially  degrade  the  mechanical  properties  of  the  glass/ceramic  sealant.  Upon  reheat¬ 
ing  to  the  SOFC  operating  temperature,  possible  self-healing  phenomenon  may  occur  in  the  glass/ceramic 
sealant  which  can  potentially  restore  some  of  its  mechanical  properties.  A  phenomenological  model  is 
developed  to  model  the  temperature  dependent  Young’s  modulus  of  glass/ceramic  seal  considering  the 
combined  effects  of  aging,  micro-voids,  and  possible  self-healing.  An  aging  time-dependent  crystalline 
content  model  is  first  developed  to  describe  the  increase  of  the  crystalline  content  due  to  the  continuing 
devitrification  under  high  operating  temperature.  A  continuum  damage  mechanics  (CDM)  model  is  then 
adapted  to  model  the  effects  of  both  cooling  induced  micro-voids  and  reheating  induced  self-healing.  This 
model  is  applied  to  model  the  glass-ceramic  G18,  a  candidate  SOFC  seal  material  previously  developed 
at  PNNL.  Experimentally  determined  temperature-dependent  Young’s  modulus  is  used  to  validate  the 
model  predictions. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  high-efficiency  energy  conversion  devices  that  are 
environmental  friendly  with  little  or  no  toxic  emissions.  The  solid 
oxide  fuel  cell  continues  to  show  great  promise  as  a  future  power 
source,  with  potential  applications  in  stationary  power  generation 
and  as  auxiliary  power  units.  Among  various  SOFC  designs,  anode- 
supported  planar  cells  have  shown  great  potential  in  delivering 
high  performance  at  reasonable  costs  [1,2].  Planar  SOFCs  offer  a 
significant  advantage  of  a  compact  design  along  with  higher  power 
densities.  In  the  meantime,  they  require  the  incorporation  of  her¬ 
metic  gas  seals  for  efficient  and  effective  channeling  of  fuel  and 
oxygen. 

Seals  are  the  most  critical  components  in  commercializing  the 
planar  SOFC  technology  [3-5].  They  must  adequately  prevent  the 
leakage  of  air  and  fuel,  effectively  isolate  the  fuel  from  the  oxi¬ 
dant,  and  insulate  the  cell  from  short  circuit.  Essentially,  there 
are  two  standard  methods  of  sealing:  compressive  sealing  and 
rigid  bonding  [6-8].  In  compressive  sealing,  a  compliant  high- 
temperature  material  is  captured  between  the  two  sealing  surfaces 
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and  compressed,  using  a  load  frame  external  to  the  stack,  to  deliver 
hermetic  sealing.  The  sealing  surfaces  can  slide  with  respect  to 
one  another  without  disrupting  the  hermeticity  of  the  seal.  This 
technology,  however,  remains  incomplete  due  to  the  lack  of  a  reli¬ 
able  high-temperature  sealing  material  [9]  and  the  difficulty  of 
designing  the  appropriate  load  frame  under  high  operating  tem¬ 
perature.  Rigid  seals  rely  on  effective  bonding  of  the  seal  material 
to  the  sealing  surfaces.  They  offer  significant  advantages  over  com¬ 
pressive  seals  which  suffer  from  problems  of  oxide  scaling  and 
chemical  stability  under  highly  reactive  environments  in  addi¬ 
tion  to  the  disadvantages  of  incorporating  an  externally  applied 
load  [6]. 

As  rigid  seal,  glass  joining  provides  a  cost  effective  and  rela¬ 
tively  simple  method  of  bonding  ceramic  and  metal  parts.  However, 
the  softening  point  of  the  glass  component  typically  limits  the 
maximum  operating  temperature  to  which  the  joint  may  be 
exposed.  As  discussed  by  Weil  et  al.  [10],  there  are  a  number 
of  other  key  materials  and  processing  variables  that  can  influ¬ 
ence  the  performance  of  glass  seals,  including  the  composition 
of  the  metal  substrate  against  which  the  seal  is  made,  operat¬ 
ing  parameters  such  as  the  expected  lifetime  of  the  device  (and 
therefore  the  seals),  as  well  as  the  degree  of  thermal  cycling 
to  which  the  seals  will  be  exposed  during  system  operation.  In 
addition,  glass  seal  in  room  temperature  is  brittle,  non-yielding, 
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and  particularly  susceptible  to  fracture  when  exposed  to  ten¬ 
sile  stresses.  Because  of  these,  much  of  the  SOFC  seal  materials 
development  effort  has  been  focused  on  developing  materials 
that  have  compatible  temperature-dependent  CTEs  for  each  of 
the  components  being  jointed,  i.e.,  the  ceramic  cell,  the  seal, 
and  the  metal  separator,  to  minimize  the  build  up  of  residual 
stresses  within  the  joint.  Only  a  handful  of  high-temperature 
glass  compositions  in  the  borate-  or  phosphate-doped  alumi¬ 
nosilicate  families  satisfy  these  requirements  [11].  The  glass  seal 
composition  used  in  this  study,  designated  as  G18  [12],  was 
developed  by  PNNL  for  planar-type  SOFC  applications.  G18  is 
a  barium-calcium-aluminosilicate  (BCAS)-based  glass  with  the 
addition  of  boron  oxide.  Silicon-based  glasses  provide  a  better 
combination  of  chemical  compatibility  and  stability  properties 
than  phosphate-  or  borate-based  glasses,  although  this  material 
is  susceptible  to  chromium  migration  when  used  with  ferritic 
stainless  steels  [13,14].  For  rigid  bond  as  such,  stress  levels  in 
the  glass  seal  and  mechanical  properties  such  as  Young’s  modu¬ 
lus  and  interfacial  strength  between  different  interfaces  become 
critical  for  long-term  reliability  of  the  seal  and  therefore  the 
stacks. 

To  obtain  a  reliable  SOFC  design  in  the  complex  operating  envi¬ 
ronments,  the  stress  level  in  the  glass  seal  materials  as  well  as 
at  the  various  interfaces  must  be  carefully  examined  and  man¬ 
aged.  To  accurately  predict  these  stresses,  accurate  material  models 
considering  the  various  aspects  of  the  seal  materials  under  dif¬ 
ferent  operating  conditions  including  aging  and  cooling  become 
critical. 

In  this  paper,  a  phenomenological  model  based  on  mechani¬ 
cal  analogs  is  developed  to  describe  the  temperature-dependent 
Young’s  modulus  of  glass  ceramic  seal  materials.  The  effects  of 
continuing  aging  and  devitrification  on  the  ceramic  phase  volume 
fraction  and  the  resulted  mechanical  properties  of  glass  ceramic 
seal  material  are  considered.  The  effects  of  micro-voids  and  self- 
healing  are  also  considered  using  a  continuum  damage  mechanics 
(CDM)  model.  The  formulation  is  for  glass/ceramic  seal  in  general, 
and  it  can  be  further  developed  to  account  for  effects  of  various 
processing  parameters.  This  model  was  applied  to  G18,  and  the 
temperature-dependent  experimental  measurements  were  used 
to  calibrate  the  modeling  parameters  and  to  validate  the  model 
prediction. 


2.  Characteristics  of  glass-ceramic  seal  materials 

The  glass  ceramic  sealant  material  examined  here,  G18,  is  a  bar¬ 
ium  calcium  aluminosilicate-based  glass  originally  melted  from  the 
following  mixture  of  oxides  (by  weight  percent):  56.4%  BaO,  22.1% 
Si02, 5.4%  A1203, 8.8%  CaO,  and  7.3%  B203  [12].  The  G18  powder  was 
milled  to  an  average  particle  size  of  20  p,m  and  mixed  with  a  pro¬ 
prietary  binder  system  to  form  a  paste  that  could  be  dispensed  onto 
the  substrate  surfaces  at  a  uniform  rate  of  0.075  g/linear-cm  using 
an  automated  syringe  dispenser.  The  glass  paste  may  be  dispensed 
onto  the  YSZ  side  of  the  bilayer  discs.  Each  disk  was  then  concentri¬ 
cally  positioned  on  a  washer  specimen,  loaded  with  a  100  g  weight, 
and  heated  in  air  under  the  following  sealing  schedule:  heat  from 
room  temperature  to  850  °C  at  10°Cmin-1,  hold  at  850  °C  for  1  h, 
cool  to  750  °C  at  5  °C  min-1 ,  hold  at  750  °C  for  4  h,  and  cool  to  room 
temperature  at  5  °C  min-1  [10]. 

Silicon-based  glasses  provide  a  better  combination  of  chem¬ 
ical  compatibility  and  stability  properties  than  phosphate- 
or  borate-based  glasses.  The  glass  wets  the  surfaces  to  be 
bonded  at  the  joining  temperature,  and  subsequent  crystalliza¬ 
tion  then  provides  mechanical  stiffness  and  strength  at  operating 
temperatures. 


Fig.  1.  Quantity  of  crystalline  phase  formed  in  glass  ceramic  sealant  (G18)  as  a 
function  of  time  held  in  air  at  750  °C:  heating  from  room  temperature  to  850  °C 
at  10  °C  min-1 ,  hold  at  850  °C  for  1  h,  cool  to  750  °C  at  5  °C  min-1 ,  hold  at  750  °C. 


2.1.  Aging  time-dependent  volume  fraction  evolution  of 
crystalline  phase 

Aging  changes  the  volume  fraction  of  the  ceramic  phase  in  the 
glass  ceramic  microstructure,  hence  influencing  seal  mechanical 
properties  such  as  Young’s  modulus,  strength,  fracture  toughness 
and  durability  [15-19].  For  glass/ceramic  sealing  materials,  the  vol¬ 
ume  fraction  of  the  ceramic  crystalline  is  aging  time-dependent 
under  the  SOFC  operating  temperature.  At  the  end  of  the  initial 
sealing  process  at  850  °C,  some  volume  fraction  of  the  ceramic  crys¬ 
talline  is  formed.  Subjected  to  the  typical  operating  environment 
of  750  °C,  the  crystallization  process  slows  down  but  does  not  stop. 
This  continuing  devitrification  causes  the  volume  fraction  of  the 
ceramic  crystalline  in  the  sealant  material  to  increase  with  the  hold¬ 
ing  time  at  the  operating  temperature.  Fig.  1  shows  the  volume 
fraction  of  various  crystalline  phases  formed  within  the  devitri- 
fying  G18  as  a  function  of  time  during  1000  h  of  aging  at  750  °C 
[10].  This  graph  was  constructed  from  a  series  of  quantitative  X-ray 
diffraction  (XRD)  measurements  conducted  on  aged  G18  samples. 

After  1  h  of  initial  sealing  process  at  850  °C,  the  crystalline  con¬ 
tents  of  BaSi03,  (Ba1.5Cao.5)Si04,  and  hexa-BaAi2Si208  are  observed 
in  the  glass/ceramic  seal.  BaSi03  is  the  main  crystalline  phase 
composition,  taking  up  approximately  40%  of  the  volume  frac¬ 
tion.  After  the  first  20  h  under  SOFC  working  temperature  of 
750  °C,  the  volume  fractions  of  both  BaSi03  and  (Bai.5Ca0.5)SiO4 
remain  relatively  unchanged,  and  the  volume  fraction  of  hexa- 
BaAi2Si208  doubled  from  7  to  15%.  With  longer  working  time, 
the  content  of  (Bai.5Ca0.5)SiO4  and  the  hexa-BaAi2Si208  began  to 
decrease,  and  the  content  of  BaSi03  started  to  increase.  A  new 
phase  of  mono-BaAi2Si208  starts  to  appear  at  about  150  h  of  work¬ 
ing  time  and  continues  to  increase  because  (Bali5Ca0.5)SiO4  and 
hexa-BaAi2Si208  are  both  transformed  into  this  phase.  The  overall 
crystalline  phase  increases  almost  linearly  with  log  (time(h))  in  the 
first  300  h  of  aging.  After  that,  the  crystallization  process  gradually 
stops  and  the  crystalline  content  in  the  material  remains  relatively 
unchanged. 

Assuming  the  devitrification  of  glass/ceramic  can  reach  an 
asymptotic  level  after  300  h  of  aging/working  time  at  750  °C,  the 
evolution  of  the  volume  fraction  of  the  crystalline  phases  with 
respect  of  time  t  may  be  expressed  as 

/c(t)=/c“-C/coo-/c0)e(_t/At)  (D 
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Fig.  2.  Comparison  of  Eq.  ( 1 )  and  experimentally  measured  crystalline  phase  volume 
fraction  for  G18. 

where  fg°  and  /°  represent  the  volume  fractions  of  the  ceramic 
phase  in  the  stable  stage  and  the  initial  stage  of  the  sealing  process, 
respectively.  At  represents  the  characteristic  aging  time.  Note  that 
various  crystalline  phases  are  lumped  into  one  term ,/c,  in  this  work. 

According  to  the  experimental  data  shown  in  Fig.  1 ,  the  param¬ 
eters  in  Eq.  (1 )  can  be  determined  for  G18  as 

/c°°  =  72.0,  /q  =  54.0,  At=  40.0(h)  (2) 

These  are  material-specific  parameters  that  are  closely  dependent 
on  the  glass  composition  and  associated  heat  treatment  processes. 
Fig.  2  compares  the  crystalline  volume  fraction  evolution  expressed 
by  Eq.  (1)  with  the  experimental  measurements  shown  in  Fig.  1. 
Clearly,  the  overall  evolution  trend  of  the  crystalline  phase  can  be 
captured  well  by  Eq.  (1 ). 

2.2.  Aging-induced  micro-voids  in  glass-ceramic 

Aging  at  high  temperature  often  leads  to  material  property 
degradations  [20,21].  One  possible  degradation  mechanism  is  the 
appearance  of  the  voids  and  cracks  [22].  For  example,  Fig.  3  shows 
the  microstructure  of  the  glass  ceramic  sealant  after  the  glass 
paste  is  sintered.  At  the  end  of  the  initial  sealing  process,  dis¬ 
tinct  boundaries  between  the  glass  (amorphous)  phase  and  ceramic 
(crystalline)  phase  can  be  observed,  and  the  fibrous  and  needle¬ 
like  crystalline  structures  can  be  seen  clearly  in  the  amorphous 
glass  phase.  Some  small  voids  and  flaws  are  also  illustrated  in  the 
microstructures. 

As  further  devitrification  takes  place  and  the  volume  fraction 
of  the  crystalline  phases  increases,  morphology  of  the  crystalline 
phases  also  changes,  leading  to  less  fibrous  needles  and  more 
diffused  phase  boundaries.  Fig.  4  shows  the  room  temperature 
microstructure  of  G18  after  aging  at  750  °C  for  1000  h.  Many  more 
small  voids  are  now  present,  particularly  in  the  amorphous  phase, 
and  the  overall  void  volume  fraction  is  also  increased  from  Fig.  3. 
The  micro-voids  are  possibly  caused  by  the  CTE  difference  between 
the  ceramic  and  the  glass  phase  during  cooling.  Different  shrink¬ 
ages  of  the  two  phases  will  likely  induce  some  voids  in  the  glass 
phase  which  has  lower  thermal  expansion  coefficients  compared 
to  the  ceramic  phase  [23-25].  Compared  with  Fig.  3,  the  distinct 
fibrous  and  needle-like  crystalline  in  the  amorphous  phase  disap¬ 
peared,  and  smeared/diffused  phase  boundaries  between  the  glass 
phase  and  ceramic  phase  are  observed  after  1000  h  aging. 


Fig.  3.  SEM  image  of  microstructure  of  G18  after  sintering  process— heat  from 
room  temperature  to  850  °C  at  10°Cmin_1,  hold  at  850  °C  for  1  h,  cool  to  750  °C 
at  5  °C  min-1 ,  hold  at  750  °C  for  4  h,  and  cool  to  room  temperature  at  5  °C  min-1 :  (a) 
magnification:  2000 x  and  (b)  magnification:  5000 x. 

To  capture  the  effect  of  aging  and  cooling  induced  microstruc¬ 
ture  changes  on  the  mechanical  property  of  the  glass  ceramic 
sealant,  a  continuum  damage  mechanics  model  [26-29]  is  used 
here  as  a  phenomenological  approach  to  the  constitutive  model¬ 
ing  of  glass/ceramic  seal.  The  microstructure-level  heterogeneous 
G18  is  represented  by  an  equivalent  homogeneous  material  with 
effective  properties.  The  model  accounts  for  the  material  damage 
because  of  various  mechanisms  (i.e.,  void  growth,  void  nucle- 
ation  and  coalescence,  decohesion  between  different  phases, 
etc.)  in  a  phenomenological  way  through  a  scalar  damage  vari¬ 
able,  D,  that  governs  the  reduction  of  the  homogenized  elastic 
modulus  [29]: 

ED(t,  T,  D)  =  E(t,  T)(1  -  D)  (3) 

where  ED  and  E  represent  the  Young’s  modulus  with  and  with- 
out  damage,  t  is  time  and  T  is  temperature.  D  is  the  quantitative 
measure  of  the  void  volume  fraction  in  comparison  with  the  virgin 
microstructure. 

2.3.  Possible  self-healing  of  glass -ceramic  seals 

Even  though  G18  exhibits  creep/flow  behavior  under  operating 
temperature,  its  behavior  at  room  temperature  can  be  character¬ 
ized  as  being  brittle.  During  thermal  cycling,  the  tensile  residual 
stresses  caused  by  the  mismatch  of  CTE  of  various  cell  components 
can  potentially  generate  many  small  cracks  and  voids  in  the  brit¬ 
tle  seal  material  after  cool  down.  If  left  untreated/unhealed,  these 
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Fig.  4.  SEM  image  of  microstructure  of  G18  after  1000  h  aging:  (a)  magnification: 
2000 x  and  (b)  magnification:  5000 x. 


small  cracks  in  the  seal  could  be  detrimental  to  the  overall  stack 
reliability  for  the  next  operating  cycle. 

In  order  to  examine  the  seal  microstructures  upon  reheat¬ 
ing,  room  temperature  Vickers  indentation  tests  and  subsequent 
reheating  is  performed  on  G18.  Fig.  5  shows  a  typical  room  temper¬ 
ature  Vickers  indentation  impression  on  the  surface  of  the  short 
crystallized  G18.  With  a  2  kg  indentation  load,  four  perfectly  sym¬ 
metrical  radial  cracks  are  generated  emanating  from  four  corners 


Fig.  6.  SEM  pictures  at  the  crack  tip  on  polished  G18  short  crystallized  with  higher 
magnification. 


of  the  indentor  edge.  Under  higher  magnification  in  Fig.  6  for  one 
of  the  radial  cracks,  it  is  observed  that  crack  propagation  followed 
a  less  tortuous  path  as  compared  to  typical  structural  ceramics 
[30,31  ].  The  white  crystallite  in  rod  shape  (likely  BaSi03 )  appears  to 
be  stronger  than  the  featureless  glass  (gray  area)  since  some  crack 
deflection  and  bridging  was  observed  as  compared  to  the  straight 
crack  path  through  the  gray  area  (likely  the  un-crystallized  resid¬ 
ual  glasses).  Nevertheless,  the  G18  showed  typical  brittleness  with 
single  crack  at  room  temperature. 

After  indentation,  the  sample  was  held  at  750  °C  for  half  an 
hour.  Figs.  7  and  8  show  the  typical  SEM  image  of  the  original 
indentation  impression  of  the  example  at  low  and  high  mag¬ 
nifications,  respectively.  It  may  be  seen  that  all  the  four  radial 
cracks  completely  disappeared,  even  under  higher  magnification. 
These  experimental  observations  can  serve  as  evidence  that  the 
glass/ceramic  sealant  material  G18  does  exhibit  some  degree  of 
self-healing:  when  reheated  again  to  high  operating  temperature, 
the  damage  in  the  glass  ceramic  sealant  may  disappear  due  to  the 
flow  of  the  glass  phase,  potentially  restoring  its  mechanical  prop¬ 
erty  to  its  undamaged  level.  Other  possible  physical  explanations 
for  this  healing  behavior  may  be  the  capillary  force  of  the  residual 
glass  phase  as  well  as  the  residual  stresses  caused  by  the  Vickers 
indentation. 


Fig.  5.  Typical  Vickers  impression  at  2  kg  indentation  load  at  room  temperature  on 
polished  G18  short  crystallized. 


Fig.  7.  SEM  pictures  at  the  impression  on  polished  G18  short  crystallized:  after 
indentation,  the  sample  was  held  at  750  °C  for  half  an  hour. 
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Fig.  8.  SEM  pictures  at  the  impression  tip  on  polished  G18  short  crystallized:  after 
indentation,  the  sample  was  held  at  750  °C  for  half  an  hour. 


3.  Assembly  of  the  phenomenological  model 

Various  formulation  and  framework  have  been  developed  to 
model  the  mechanical  behavior  of  composite  materials  [33-38].  For 
simplicity,  we  adopt  here  the  simple  Rule  of  Mixture  approach.  In 
general,  the  Young’s  modulus  of  the  multi-phase  composite  mate¬ 
rials  can  be  expressed  as  [39]: 

m 

*  =  5>  (7) 

1=1 

where//  and  Ez  represent  the  volume  fraction  and  modulus  of  the  ith 
phase,  and  m  is  the  number  of  the  phases  in  the  composite  material. 
The  volume  fractions  of  all  the  phases  should  satisfy: 

m 

X/  =  1  (8) 

i=l 

In  the  case  of  the  two-phase  glass/ceramic  sealant  composite  mate¬ 
rial,  the  above  equation  can  be  rewritten  as 


It  should  be  mentioned  that  the  concept  of  self-healing  glass 
seals  was  first  reported  by  Singh  [32].  The  rationale  behind  this 
concept  is  that  at  the  SOFC  operating  temperature,  a  sealing  glass 
with  appropriate  properties  can  heal  cracks  created  during  thermal 
transients.  The  advantage  of  this  approach  is  that  materials  with 
dramatically  different  expansion  can  potentially  be  used  for  seals, 
because  at  the  cell  operating  temperatures,  CTE  mismatch-induced 
thermo-mechanical  stresses  can  be  relaxed  out.  If  proven  feasible, 
this  type  of  seals  can  tolerate  some  of  the  thermal  expansion  mis¬ 
match  and  still  form  a  functioning  seal  between  materials  with 
significant  expansion  mismatch.  Flowever,  there  are  still  a  num¬ 
ber  of  challenges  in  making  a  functioning  seal  with  self-healing 
glasses.  The  first  challenge  is  that  the  glass  seal  must  maintain  con¬ 
tact  with  all  the  components  to  be  sealed  without  excess  creep/flow 
at  cell  operating  temperature  of  750  °C.  Concepts  involving  various 
ceramic  phase  stoppers  have  been  proposed  to  address  this.  In  addi¬ 
tion,  the  self-healing  glass  seal  must  be  able  to  remain  chemically 
stable  under  SOFC  operating  temperature  for  the  designed  stack 
operating  time. 

For  SOFC  stacks  under  deep  thermal  cycles,  cracking  and  micro¬ 
voids  may  occur  in  the  glass  ceramic  sealant  due  to  the  tensile 
residual  stress  caused  by  CTE  mismatch  of  different  components 
during  cooling.  Flowever,  upon  reheating,  the  damage/crack  may 
be  healed  at  temperatures  above  some  threshold  temperature  due 
to  the  flow  behavior  of  the  glass  phase.  Therefore,  the  scalar  dam¬ 
age  D  defined  in  Eq.  (3)  is  temperature  and  aging  time-dependent, 
which  may  be  simplified  as  the  product  of  the  following  uncoupled 
factors: 


D=AdDt  (4) 

Here  AD  is  the  aging  influence  function,  representing  the  increase 
of  micro- void  volume  fraction  with  aging  time: 

Ad  =  1  -  e-fo  (5) 


£=/c£c  +  (l -/c)£g  (9) 

or 


E=fc(Ec-EG)  +  EG 


(10) 


where  subscripts  “C”  and  “G”  refer  to  glass  phase  and  ceramic 
phase,  respectively.  In  general,  the  modulus  of  the  glass  phase  is 
temperature  dependent,  and  modulus  of  the  ceramic  crystalline  is 
independent  of  temperature. 

Substituting  the  aging  time  dependent  volume  fraction  of  the 
ceramic  phase,  i.e.,  Eq.  (1),  into  Eq.  (10)  leads  to: 

E(t,  T )  =  t /c°°  -fgJ-t/At)XEc  -EC(T))  +  EC(T)  (11) 

Incorporating  the  effects  of  micro-damage  and  self-healing 
described  in  Eqs.  (3)-(6)  into  Eq.  (11)  yields  the  following  phe¬ 
nomenological  model  describing  the  elastic  modulus  of  the  glass 
ceramic  materials: 


E°(t,  T,  D)  =  [(fc°°  -/Oet-^XEc  -  EC(T))  +  EC(T)] 


X 


1  -  ^(1  -  e~t/tc 


n 

2.0  -  arctan 


(12) 


Again  /°°  and  /°  represent  the  volume  fractions  of  the  ceramic 
phase  in  the  stable  stage  and  at  the  initial  stage  of  the  sealing  pro¬ 
cess,  respectively.  At  represents  the  characteristic  aging  time.  D0  is 
the  damage  parameter  at  room  temperature,  and  Tth  refers  to  the 
threshold  temperature  above  which  damage  of  the  glass  ceramic 
composite  starts  to  heal.  R  is  a  material  constant,  t  is  the  aging 
time,  and  tc  represents  the  characteristic  time  for  aging  induced 
micro-damage. 


4.  Comparison  between  experimental  measurements  and 
predictions 


where  t  is  the  aging  time,  and  tc  represents  the  characteristic  time. 

Dj  in  Eq.  (4)  represents  the  temperature-dependence  of  dam¬ 
age,  reflecting  the  possible  self-healing  behavior  of  G18  above  the 
threshold  temperature,  Tth: 

D(T)  =  D0  */2-0-arctan[(r-rth)/R]  (6) 

7 T 

where  D0  is  damage  parameter  at  room  temperature,  and  Tth  is  the 
threshold  temperature  above  which  damage  of  the  glass  ceramic 
composite  starts  to  heal,  and  R  is  a  constant. 


4.1.  Temperature-dependent  modulus  measurements 

The  dynamic  resonance  technique  (ASTM  C1198)  was  used  to 
measure  the  elastic  moduli  of  the  specimens  of  both  non-aged 
and  1000  h  aged  short  crystallized  G18  at  room  and  elevated 
temperatures  [40].  The  instrument  used  to  measure  the  resonant 
frequencies  was  constructed  following  ASTM  C1198.  The  specimen 
was  suspended  within  a  resistance-heated  furnace  using  single¬ 
crystal  sapphire  fibers  with  beaded  ends.  Small,  narrow  notches 
were  cut  into  the  corners  of  the  specimen  to  accept  the  fibers. 
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Table  1 

Young’s  modulus  aged  and  un-aged  G18  at  room  and  elevated  temperatures 


The  fibers  coupled  the  specimen  to  two  piezoelectric  transduc¬ 
ers  mounted  above  the  furnace.  One  transducer  was  used  as  a 
driver,  and  the  other  as  a  receiver.  A  computer-controlled  sys¬ 
tem  was  used  to  send  frequency  signals  sweeping  from  20  Hz  to 
20  kHz  through  the  driver  transducer.  The  signal  from  the  receiv¬ 
ing  transducer  was  fed  back  into  the  computer,  and  the  resonant 
frequencies  were  identified  and  recorded.  Both  the  transverse 
and  longitudinal  resonant  frequencies  could  be  reliably  detected 
at  the  corners  of  the  specimen.  Resonant  frequency  data  were 
then  used  in  calculating  the  elastic  moduli  with  the  following 
equations  [41]: 


E  = 

G  = 


mf }  L3 

0.9465 

b  h3 
4Lmft  B 


bh  1  +A 


(13) 

(14) 


where  E  and  G  are  the  Young’s  modulus  and  shear  modulus,  respec¬ 
tively;  m  is  the  mass  of  the  bar,  b ,  h  and  L  are  the  width,  thickness, 
and  length  of  the  bar,  respectively.  /f  represents  the  fundamental 
resonant  frequency  of  the  bar  in  flexure;  and  T\  denotes  the  correc¬ 
tion  factor  for  fundamental  flexural  mode  to  account  for  the  finite 
thickness  of  the  bar  and  Poisson’s  ratio,  etc.  The  Young’s  and  shear 
moduli  were  then  used  to  calculate  Poisson’s  ratio: 


M  = 


(15) 


Temperature  (C) 

Fig.  9.  Young’s  and  shear  moduli  for  aged  and  un-aged  G18  at  room  and  elevated 
temperature:  the  furnace  was  heated  at  a  rate  of  approximately  12°Cmin_1,  and 
was  held  for  15  min  at  each  temperature  for  equilibration. 


400  to  600  °C.  Above  600  °C,  the  moduli  start  to  decrease  with 
increasing  temperature. 

For  temperatures  lower  than  400  °C,  the  measured  modulus  for 
the  aged  glass  ceramic  sealant  is  consistently  less  than  that  of  the 
non-aged  (short-term  aged)  glass  ceramic.  At  about  500 °C,  the 
moduli  for  the  aged  and  non-aged  samples  cross  over.  For  temper¬ 
atures  above  600  °C,  the  moduli  trend  reverses  itself:  the  modulus 
for  1000  h  aged  G18  is  higher  than  that  of  the  non-aged  G18. 

This  distinct  temperature-dependent  modulus  behavior  for  the 
aged  G18  is  another  motivation  for  our  current  work.  Our  goal  is  to 
see  whether  this  unique  behavior  can  be  predicted  by  coupling  the 
various  phenomena  such  as  aging,  cooling  and  reheating  induced 
self-healing.  The  eventual  goal  is  to  develop  the  predictive  capa¬ 
bility  for  glass  ceramic  seals  such  that  better  seal  materials  can  be 
developed  for  improved  reliability  and  durability. 


The  precisions  of  the  E ,  G  and  /z  values  were  calculated  by  propa¬ 
gating  the  uncertainty  of  each  term  through  the  equations. 

The  test  specimen  was  a  right  parallelepiped  measur¬ 
ing  50.75  mm  x  25.34  mm  x  5.96  mm.  All  surfaces  were  ground 
smooth.  Resonant  frequencies  were  measured  at  room  temperature 
and  in  100  °C  intervals  from  200  to  800  °C.  The  furnace  was  heated 
at  a  rate  of  approximately  12  °Cmin-1,  and  was  held  for  15  min  at 
each  temperature  for  equilibration. 

The  measured  Young’s  modulus  for  the  non-aged  and  1000  h 
aged  G18  at  room  and  elevated  temperature  is  tabulated  in  Table  1. 
Both  of  the  measured  Young’s  modulus  and  shear  modulus  for 
the  non-aged  and  1000  h  aged  G18  at  room  and  elevated  tem¬ 
perature  are  shown  in  Fig.  9.  For  non-aged  short  crystallized 
glass/ceramic  sealant  (G18),  both  the  Young’s  modulus  and  shear 
modulus  decrease  with  increasing  test  temperature.  The  resonant 
peaks  became  very  broad  at  about  700  °C,  making  it  impossible  to 
differentiate  between  the  longitudinal  and  transverse  peaks  with 
any  certainty.  This  may  have  been  caused  by  the  softening  of  the 
glass  phase  present  in  these  specimens.  Additional  measurements 
taken  again  at  room  temperature  (30  °C)  after  the  high-temperature 
tests  achieve  the  same  moduli  as  the  pre-test  values. 

The  measured  Young’s  modulus  and  shear  modulus  for  the 
1000  h  aged  G18  show  some  very  distinct  and  interesting  features. 
For  temperatures  lower  than  400  °C,  both  Young’s  modulus  and 
shear  modulus  increase  with  increasing  test  temperature.  The  mod¬ 
uli  exhibit  somewhat  temperature  independent  behaviors  from 


4.2.  Modeling  results 

For  the  phenomenological  model  described  in  Eq.  (12),  the 
modulus  of  the  ceramic  phase  is  considered  to  be  temperature 
independent,  and  modulus  for  the  glass  phase  is  considered  to  be 
temperature  dependent.  In  general,  the  Young’s  modulus  of  glass 
drops  dramatically  over  the  glass  transition  temperature,  Tg  [42]. 

The  input  Young’s  moduli  used  for  the  glass  and  crystalline 
phases  are  shown  in  Fig.  10  in  dashed  lines  [42].  The  predicted 
Young’s  moduli  of  G18  under  different  aging  times  are  shown  in 
solid  lines  in  Fig.  10.  The  following  parameters  for  Eq.  (12)  are 
calibrated  based  on  the  experimental  data: 

D0  =  0.24,  Tth  =  375  °C,  R  =  50  °C,  tc  =  100h 

With  the  proposed  phenomenological  model,  reasonably  good 
comparisons  between  the  predicted  and  measured  Young’s  modu¬ 
lus  have  been  achieved  for  both  the  aged  and  un-aged  G18.  This  is 
particularly  true  for  temperatures  lower  than  600  °C.  Above  600  °C, 
the  predicted  modulus  decreases  more  rapidly  with  temperature 
than  the  measured  data.  This  is  because  the  predicted  modulus 
is  very  sensitive  to  the  input  temperature-dependent  modulus 
for  the  glass  phase.  For  temperatures  above  600  °C,  the  Young’s 
modulus  for  glass  may  vary  dramatically  [43,44].  The  exact  high- 
temperature  Young’s  modulus  for  regular  glass  is  rarely  reported, 
and  the  high-temperature  Young’s  modulus  for  the  glass  phase  in 
G18  has  not  yet  been  separately  determined.  Therefore  for  temper- 
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atures  above  600  °C,  the  glass  moduli  used  in  Fig.  10  could  have 
been  much  lower  than  the  actual  values. 

It  is  interesting  to  note  that,  by  considering  the  combined  effects 
of  aging,  micro-voids  and  self-healing,  the  distinct  feature  of  the 
moduli  ‘cross-over’  can  be  predicted  by  the  simple  phenomenolog¬ 
ical  model:  Below  500  °C,  the  1000  h  aged  sample  has  lower  Young’s 
modulus  compared  to  the  4  h  aged  sample;  Above  500  °C,  the  trend 
reverses  itself  when  considering  self-healing  of  the  glassy  phase. 

5.  Conclusions 

We  study  the  temperature-dependent  Young’s  modulus  of  G18, 
a  glass  ceramic  material  for  SOFC  sealing  applications.  The  com¬ 
bined  effects  of  aging,  cooling-induced  micro-voids  and  possible 
reheating  related  self-healing  are  examined  and  modeled  with 
some  simple,  phenomenological  models.  After  the  initial  seal¬ 
ing  process,  the  crystallization  process  of  G18  slows  down  but 
does  not  stop,  and  the  volume  fraction  of  the  crystalline  phase 
increases  with  aging/working  time.  In  addition,  aging  causes  the 
diffusion/smearing  of  the  boundaries  between  the  crystalline  and 
the  amorphous  phases  and  potentially  changes  the  mechanical 
properties  of  the  amorphous  phase. 

Upon  cooling  to  room  temperature,  shrinkage  micro-voids  will 
be  formed  because  of  the  CTE  differences  between  the  crystalline 
and  the  amorphous  phases.  On  the  stack  level,  cooling  induced 
micro-cracks  are  also  likely  to  occur  due  to  the  CTE  mismatch 
between  different  stack  components.  These  micro-voids  and  micro¬ 
cracks  can  noticeably  degrade  the  Young’s  modulus  of  the  glass 
ceramic  seal  at  room  temperature.  However,  when  reheated  back  to 
SOFC  working  temperature,  this  study  shows  that  G18  does  exhibit 
some  degree  of  self-healing  behaviors  because  of  the  flow  char¬ 
acteristics  of  the  glass  phase  at  high  temperatures.  Therefore  the 
mechanical  property,  i.e.,  Young’s  modulus  of  the  glass/ceramic  seal 
material  can  be  potentially  restored  to  its  undamaged  level  at  high 
temperatures. 

The  phenomenological  model  based  on  some  simple  mechani¬ 
cal  analogs  is  developed  to  capture  the  above-described  mechanical 
behaviors  of  glass  ceramic  seal  materials.  The  aging  time  depen¬ 
dent  crystalline  phase  evolution  model  was  first  developed  to 
describe  the  increase  of  crystalline  content  due  to  the  continuing 
devitrification  during  operation.  A  continuum  damage  mechanics 
model  was  adapted  to  model  the  effects  of  micro-voids  and  self- 
healing.  Reasonably  good  comparisons  between  the  measured  and 
the  predicted  temperature-dependent  Young’s  modulus  have  been 
obtained.  The  modeling  parameters  presented  here  are  calibrated 
with  the  experimental  data  for  G18,  yet  the  modeling  framework 


should  be  applicable  to  glass/ceramic  seal  in  general.  The  model 
can  be  further  developed  to  account  for  the  effects  of  various  pro¬ 
cessing  parameters  on  the  mechanical  properties  of  glass/ceramic 
seal  materials. 
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